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Abstract 

Nonlinear solutions for the large-scale flow of interstellar gas in the presence 
of a spiral gravitational field have received considerable attention because this 
model allows to explain a lot of observations for galaxy disks. However this 
investigations were forced to limit by numerable analysis of the problem because 
of its nonlinearity. 

In this paper we wish to carry out the analytical expression which allows 
us to describe the super sonic nonlinear interstellar gas flow in rapidly rotating 
galaxy disk which is made by the spiral density wave. 

One of characteristic parameters of theory is the amplitude of spiral density 
wave potential corresponding to separatrix, which separates super sonic flows 
from flows containing the jump from super sonic flows to subsonic. 

We have defined the dependence of perturbing potential value, under which 
the galaxy shock wave appeared with respect to the parameters characterizing 
the gas disk (such as sound speed in gas, disk rotation speed, ’’spiral design” 
rotation speed). 

In recent years there have been several numerical simulations for nonlinear solu¬ 
tions of large scale flow of interstellar gas in the stellar gravitational field of spiral 
galaxies (Lin, Yuan & Shu 1969; Roberts 1969; Woodward 1975) by using the wave 
theory of spiral structure of galaxy. This concept has played a remarkable role in un¬ 
derstanding of many processes occurring in the galaxies such as the structure of dust 
lanes, observed along the inner edges of spiral arms in many galaxies (Fujimoto, 1966; 
Lynds, 1970), the enhanced synchrotron radiation from spiral arms (Mathewson et 
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al, 1972; Sofue, 1985), the radio emission of HI at the wavelength A = 21m (Berman 
& Mishurov 1980), a trigger mechanism for star formation, creating the narrow bands 
of young highly luminous stars which delineate the spiral arms (Roberts, 1969; Shu et 
ai, 1972; Shu, 1973) and etc. It was shown that the spiral density wave, propagated 
in the stellar galaxy disk and interpreted as the galaxy spiral arms, is induced large 
nonlinear perturbations in the axial symmetric gas flow of the galactic disk. Such 
nonlinear phenomena take place in the gas even though the amplitude of the spiral 
stellar held is small. This is because the response to a spiral gravitational potential 
induced by the stellar density wave is roughly proportional to a“^, where a is the 
velocity dispersion, for stars, or sound speed for the gas (Woodward 1975). For the 
gas a ~ 8 while for the stars a ~ 40 So, if the amplitude of the perturbation 
held is a few per cent and we can use the liner theory to calculate a disturbance in the 
stellar disk, the perturbation of the gas disturbance is much stronger and we must use 
nonlinear theory. Nevetherless, the nonlinearity makes this model more contextual 
in the sense of physical ideas, which let to explain the morphology peculiarities of 
building and evolution of galaxies. 

In this paper we wish to carry out the analytical expression which allows us to 
describe the disturbance of the super sonic nonlinear interstellar gas how in rapidly 
rotating galaxy disk which is made by the spiral stellar density wave. 

The hydrodynamic equations for the steady nonlinear gas how (self-similar solu¬ 
tions) in the frame rotating with angular velocity Qp are, at position r, 

V • (pv) = 0 (1) 


(v ■ V)v -|- 2Q,p X V 


--VP-V(<I>-^H/r2) 
P ^ 


( 2 ) 


where p, v and P = a^p (we have assumed that the gas is isotermal, i.e. a 
is isothermal sound speed) are respectively the density, velocity, and pressure of the 
interstellar gas; $ is the gravitational potential of the stars contained an unperturbed 
axisymmetric part as well as non axisymmetric perturbation induced by spiral wave; 
r is a distance from the centre of the galaxy disk. In the absence of a stellar density 
wave the gas how is supposed to be purely circular, so the gravitational potential <F is 
axisymmetric and can be determined by the condition <F(r) = where H = r2(r) 

is an angular velocity of rotating galaxy disk at position r 

The unperturbed axisymmetric solutions of equations (1) and (2) for rapidly ro¬ 
tating thin disk are 


VzQ 0, Vro 0, 'r(f2 f^p) 

where, Vro,v^o,'Vzo are in cylindrical polar coordinates. 
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Describing the nonlinear gas flow perturbed by the spiral stellar density wave each 
hydrodynamic variable can be written as a combination of an unperturbed axisym- 
metric part and a perturbation, denoted by subscripts 0 and 1 respectively (Roberts 
1973) then the gravitational potential is <l> = $0 + According to other authors, 
beginning from Roberts (1973) it can be suggested that the perturbing potential <h 
is similar to the spiral. It is then useful to introduce the spiral coordinates z). 
The spiral coordinates is the set of coordinates fixed in the rotating frame that are 
parallel and perpendicular to the spiral equipotential square so rj is constant along 
the spiral arms while ^ is constant everywhere on lines orthogonal to the spiral arms 
and 2 : is constant everywhere on plane parallel to galaxy disk. This coordinates are 
related to the plane cylindrical polar coordinates (r, 0) by 

2 

rj = -Inr + 20 

tani 

2 

f = —2 Inr H-0 

tani 

where i is angle of inclination of a spiral arm to the circumferential direction; this angle 
is related to the radial wavenumber of the spiral wave k by conditions tani = 
and /c = — A is the wave length, so we get 


df] = —kdr + 2d0 
dv 

di =-2 - krd(t) 

r 

The unperturbed gas velocity components parallel to the ((^, rj, z) directions re¬ 
spectively are 


uyo = (D — Dp)r cos i, v±o = (D — Dp)r sin i 



Qp) cosi, 


where ny, v± are the gas velocity components parallel and perpendicular to the 
spiral arms respectively. It is easy to see that ^ ~ A. Since each hydrodynamic 

Sr' Sr' J J /-T’ J J 

variable can be written as a combination of an unperturbed part and a pertubation 
we can define as r’y = vyo + = vj_o -l- nj_i, p = po + pi 

According to Roberts (1969) we can assume that the spirals are tightly wound, 
i.e. f 1, it corresponds to the short, comparing with scale r, waves, i.e. A r, and 
allows to use WKB approximation for the spiral wave description. We would like to 
note that the spirals in this case appear to be ’’tightly” wound. Then it can be shown 
(Nelson & Matsuda 1977), that the gradients for the spiral perturbations satisfy 


d 

TT 7 < sm 


A 

dr)^ 
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and it can be also found that for any function ^(r) with scale length r i.e. 
we get 


^ ~ sin z ^log(n||i,n_Li,pi), 

In this approximations it can be neglected derivatives of perturbed quantities in 
the direction of the ’’slow” coordinate ^ along the spiral arms with respect to their 
derivatives in the direction of the ’’rapid” coordinate r] which is perpendicular to the 
spiral arms. In addition, it can be said that the typical length in this problem is 
I ~ A, besides that j 1 (here A is a typical thickness of the galaxy disk at 
position r ~ lOkpc) than because of equation (1) Vz = 0 (it is right because in such 
conditions the last term in equation (1) is dominant one) Finally assuming the 
perturbation potential as 


^i(h) 


cos^ i 

F -;-COS f] 

k 


where F is the amplitude of the perturbation potential due to the stellar density 
spiral wave as a fraction of the unperturbed axisymmetric gravitational potential. 

It can be introduced scale = 11.34^ for ny as well as scale for vi_ = 18.1^ 

than uin—^ = —z/ = 0.722 Here are used the same set of parameters as 

used by Shu et al (1973): H = 24.7 kip = 13.5 at position r = 10 kpc ; 

X = 2fi((l + ^f cos2 0)^ = 31.0^ 

Eliminating pi between equations (1) and (2), according to Shu et al (1973) we 
can write this equations in dimensionless form 


(— 1 / + uV — c du „ . 

- — = V — f sin 71 

—v + u dr] 

, , dv 

u + {-u + u)— = 0 
dr] 


(3) 


where u, v are the dimensionless velocity variables perpendicular v± and parallel uy 

to the spiral arms respectively; / = F^^kr ~ 0.1 —0.2 is the dimensionless amplitude 
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of the stellar density spiral wave; c = ^^ = 0.195 is the square of the dimesionless 
sound speed 

Equations (3) have been solved numerically by other authors mentioned above. 
Eliminating v between these equations we can get a equation which is to be solved 
for u 


{— 1 / + uY — c d^u 
—u + u dr]‘^ 


(-Z/ + -u)^ + c 

{—u + uy 


,du.n 


u 


— p + u 


—/cos r] 


(4) 
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Beginning with some value of / the solution of equation (4) contains an irregularity 
is due behavior of the equation solution in neighborhood some point aX u = v + 

If we introduce a new variable by condition u = u — v — ^Jc then we get 


(2y/c + n) 

777—F-TT'UM 

2(-\/c + u) 


+ 


(1 


2y/c -\- U 

2(y/c + 12)2 / 


-/2 


u 


2^c{^c + uf 2 ^^ + 


^ + / cos rf) 
VC 


In the flow zone, which lies in the neighborhood of this point, gas parameters 
corresponding to the explosive solution are exposed by sharp changes, and along the 
spatial coordinate a narrow transitional area appears, in this area gas parameters 
constantly change: a so called jump of density if in the chosen frame this area is 
stable, or a shock wave if the transitional area is transferred in space with respect to 
time. Nevertheless, for the counting of flow parameters in this area we need to take 
into account the additional non-gas dynamic transfer of impulse and energy, which 
corresponds to the viscosity and heat of dissipation. The latter may be taken into 
account if we add in the simplest case ’’viscous pressure” to the gas-dynamic pressure. 

By assuming in (4) that u = fi/y and — c)~^ we get 


y” -I- u‘^y = up' cos r] -|- fiplv^uP -|- pyy” -|- {2v^uP — 1)?/'^ -|- 

-I- uPiyp' — 2ycosr])) + f'^Pu'^{—3y‘^y" — 2yy'^ + y"^ cosy) + 

+ fW{y^" + y^y'^) 


If we choose the terms in the same powers of / we get the following set of equations 
which are to be solved successively for yi; y 2 ] ys] ■ ■ ■ 


yo” + u'^yo = cos y 

yi + U)‘^yi = i2Puj^ + l)yoyo' + - 1) + ^^{yo - 2yo cos y) 

yP + u}^y 2 = {2Pu^+ l){yoyi”+ yiyo”) + { 2 ^ 0 }"^-l)2yoyi + 

+ uj‘^{2yoyi - 2yi cos y) + Puj‘^{-3ylyP - 2yo?/o'^ + yl cos y) 


yP + u^ys = {2Pu‘^ + 1) {y^yP + yiyi" + ?/ 22 /o") + 

-F (2z/^ct;^ - + 2yoy2) +uj‘^{yl + 2yo?/2 - 2?/2COS7]) -h 

-|- P{—3{yQyi + 2yo?/iW0 ~ 2(2^o?/o^?/i^ + yWo ) + 

-h 2yo?/i cos y) + Pu‘^{ylyo' + ?/o 
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This set of equations has been used by Shu,Milione & Roberts (1973) for discussing 
the slightly nonlinear regime and comparing it with the results of the linear theory by 
Lin et al. (1969). At the same time the authors noted that for such value as realistic 
gravitational potential, the series which are solutions of this equation set may converge 
very slowly, if at all. So they used a numerical method of investigation this problem 
and they got that for the typical value of parameter /, such as /c = 0.105 the behavior 
of the solution was changed and if / was being increased above fc the solution could 
not remain smooth everywhere and had to contain a jump. Physically, the jump 
corresponds to a shock wave which are formed at the transition from supersonic gas 
flow to subsonic. 

The solution y{r]) can be written in the form 

OO OO 

yiv) = Ai,k cos{kr]) (5) 

i=0 k=0 

In summation over i the dominant contribution arising from a first few terms 
because if i increases the series coefficient decreases very quickly, on the contrary, 
in summation over k dominant contribution is arisen from the remote terms of this 
series because the coefficient Ai^k is changed with respect to k very slowly. The 
necessity to take into account the remote terms is so important as the gas flow velocity 
is more closely to sound speed, this corresponds with that the value of parameter / 
is approached to its typical value mentioned above. 

The main coefficient A^ n can be obtained from solution of n — 1 equation, to do 
it we must join terms which are proportional to cos(n + l)r], if n > 0 we get 


Un 


(n+1) 

[u u Si{n) + -S 2 {n) 


u 


{n + 1)^ — 


{n + 1) 


-(S'i(n) — ttn-l) + 


,Si{n — 1) 
4 ^ (n + 1)2 




SSsin) - 2 S' 4 (n)) H- —{ 85 ( 71 ) + Se{n))) cos(n + l)r] + 


( 6 ) 


the terms contained cos((n — l)r]) and etc. 


Where 


Si{n) 


n —3 

pQ- n —1 T ^2k=0 ^kO'n—l—k 

X)a:=0 ^k^n—l—k 


if n is odd 
if n is even 
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if n is odd 
if n is even 


<52 (n) 


n-2k-l\2„ „ 

) (^kOin—l—k 

I n-2k-l \2 „ „ 
Z^A:=0v n+1 ' ^k^n-l-k 


T.kl 

n-2 


n—2n—2—k 


Ssin) = J2J2 

k=0 m=0 


n 


m 


n + 1 


fakaman-2-k-m 


04 ( 71 ) ( / , / , ^ ^^2 (^k^m(^n—2—k—m 


k=0 m=0 


n—3n—3—kn—3—k—m _ o _ L _ ™ _ / 

*55(n) = E E Y. i- 

k=0 m=0 1=0 


n + 1 


) ^k^m^l^n—3—k —m—l 


n—3n—3—kn—3—k—m 


Soin) = E E E 

k=0 m=0 1=0 


(m + l)(n — 2 — k 


m 


1) 


{n + 1)^ 


^k^m^l^n—3—k —m—l 


where is a coefficient at cos(?7, + l)r] in the solution yn which is related to Ao^n 

2 

by relation 24o,n+i = Recursion relation for this coefficients can be 

obtained by using the successive approximations method to solve the equation (6) for 
n approach. The hrst approach can be written in form 

a„^.-(1.89/ao(.V + i))-n (, P) 


For large n, taking into account that the product approximates rapidly, it can be 
substituted by limit nm=i(l — (equality u to whole number corresponds 

to the Lindblad’s resonance). So for series (5) the following expression can be written: 

, ,2 3 00 1 2 

y-^ —-cos(? 7 ) - 0.0342^-^(1.89/(i^^ct;^ +-)^—-)''cos {k + l)r] (8) 

— 1 sm TTU r~, 2 — 1 


The numerical coefficients appeared in the last two expressions are the results of aver¬ 
aging the difference between the product of hnal number of the terms ]Xin=o 
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and its limit when the number is inhnite and the result of the calculation of the sums 
S'i(n), S' 2 (n) and etc. The latter can be done while using the method of continued 
fraction. 

Trigonometric series in ( 8 ) in condition that | 1.89/(i^^ci;^ + |< 1 converges 

and, according to Prudnikov, Brychkov & Marichev (1981) can be expressed in ele¬ 
mentary functions, so we have: 


y - 



cos(? 7 ) — 0.0342 


TTU^ COS f] — a 

sin TTu 1 — 2 q ; cos r] + 


(9) 


Here 

a = 1.89/(z/2ct;2 + ]-) ^ 

Using (9) we will put down at last the expression which consists the main part of 
speed component perpendicular to wave front, i.e. taking into account that (9) 
to the large extend is adequate to large k from (5), we will change in (9) at least the 
hrst member by the exact solution of the equation for the hrst approximation, so we 
get: 


+ 4r(' 


2 - 1 


r-f 


u 


- 1 


cos 1 ] 


+/ 


— 1 — l){u^ — 4) 2 


+0.0342/ 


TTLJ 


COS 7] — a 


sin TTcj 1 — 2a cos 77 + 0 ;^ 


( 10 ) 


Integral curves (10) of equation (4) for several values of parameter / are given on 
Figure 1. As we can see from this picture with parameter / increasing, i.e. with the 
growth of perturbed potential amplitude the gas flow speed, at 77 = 0 , is decreasing 
and may reach sound speed in gas. The integral curve, corresponding to this condition, 
has a sharp jump if 77 = 0 and presents a separatrix which separates the solutions 
corresponding to the gas flow with supersonic speed from the solution, in which there 
are areas of both supersonic and subsonic gas flow. If one prefers to see a results 
represented as the velocity plane its pictured at Figure 2. 


The amplitude value of perturbed potential fc which corresponds to this separatrix 
was dehned by Shu et al. numerable. Formula (10) allows to dehne the dependence fc 
respect to the parameters of the rotating gas disk. For this we take in (10) u± = \fc 
if 77 = 0 , then solving the equation in respect to / we get: 


fc 


1 sin TTo; ^ 1 

2 0.0342 -0.0616(1 - 1.89(1 + ^)(z/2o^2 + |)) 






0.12 



Figure 1. The phase portrait of equation. Curve A assumes the amplitude of the 
perturbation potential such as / = 0.11; B assumes / = 0.10; C assumes / = 0.095; 
D assumes / = 0.085; E assumes / = 0.055; The sound speed is marked by broken 





Figure 2. The supersonic flow in the velocity plane. Curve A assumes the amplitude 
of the perturbation potential such as / = 0.11; B assumes / = 0.10; C assumes 
/ = 0.095; D assumes / = 0.085; E assumes / = 0.055; For / = 0.11 the flow 
contains a cusp at the location of the minimum of the spiral potential, r; = 0. The 
unperturbed axisymmetric velocity is marked by a cross. 
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Here the same parameters with which Shu et al (1973) solve equation (3) numerable, 
and the value we get is practically coincides with their result. 

The author thanks Professor Yu.N.Efremov & Professor A.D.Chernin from Stern¬ 
berg Astronomical Institute, Moscow State University for suggesting this problem as 
well as for their interest to this work. The author also thanks Professor Yu.N.Mishurov 
from Rostov State University, Rostov-on-Don for giving the useful references. 
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